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Excitation functions for total reaction cross sections, gr, were measured for the light, mainly 
proton-rich nuclei 6 Li, 7 Be, 10 B, 9 ' 10 ' U C, 12 N, 13 ' 15 0, and 17 Ne incident on a Si telescope at ener- 
gies between 15 and 53 MeV/nucleon. The telescope served as target, energy degrader and detector. 
Proton-removal cross sections, <T2 P for 17 Ne and a p for most of the other projectiles, were also mea- 
sured. The strong absorption model reproduces the ^-dependence of <tr, but not the detailed 
structure. Glauber multiple scattering theory and the JLM folding model provided improved de- 
scriptions of the measured or values, rms radii, extracted from the measured or using the optical 
limit of Glauber theory, are in good agreement with those obtained from high energy data. One- 
proton removal reactions are described using an extended Glauber model, incorporating second 
order noneikonal corrections, realistic single particle densities, and spectroscopic factors from shell 
model calculations. 

PACS numbers: 24.10.-i,25.60.-t, 25.60.Dz, 25.60.Gc, 25.70.Mn 



I. INTRODUCTION distributions of nuclear matter, shell structure, excited 

states, and decay modes. Often, the one-particle sepa- 
Thc interesting properties of radioactive nuclei pro- ration energies are small, of the order of 1 MeV or less, 
duced in the laboratory include their lifetimes, sizes and 
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These small separation energies lead to a wealth of phe- 
nomena including soft collective modes, exotic transition 
strengths between low-lying states, changes in shell struc- 
ture, long-tailed density distributions, and, perhaps most 
dramatically, halo nuclei. In onc-particlc removal reac- 
tions, the parallel and transverse momentum distribu- 
tions of the core-like fragments are very narrow as com- 
pared with those of normal nuclei, indicating increased 
nuclear size. The total reaction and breakup cross sec- 
tions are large, also reflecting the increased nuclear size. 
Spontaneous 2p radioactivity has been observed recently 
near the proton drip line0,|3|- The increased density in 
the tails of the matter distribution generates a compe- 
tition between the increased refractive power of the real 
optical potential and the increased absorption due to the 
imaginary part, leading to exotic shapes in heavy ion 
elastic scattering angular distributions. 

Total reaction cross sections, or, and proton-removal 
cross sections, cr p , contain complementary information 
about the size and matter distributions of atomic nuclei. 
Various reactions at all impact parameters contribute to 
or, which therefore reflects mainly the rms nuclear ra- 
dius. Breakup, however, is a peripheral process, hence 
a p is sensitive mainly to the surface distribution. There- 
fore these combined data are needed to better describe 
the matter distribution, and the best test of theory is 
obtained by fitting them simultaneously. 

Measurements of both or and o p can be used to iden- 
tify proton-halo nuclei. Among the light nuclei, 8 B has 
been identified as a proton-halo nucleus. The evidence 
includes its enhanced or at low energies , as well as a 
narrow longitudinal momentum distribution of the 7 Be 
core following proton removal. [3- 0- 

We report or measurements for the stable and short- 
lived light nuclei 6 Li, 7 Be, 10 B, 9 ' 10 ' U C, 12 N, 13 < 15 0, and 
17 Ne incident on Si targets at energies ranging from 15 
to 53 MeV/nucleon. Most of these nuclei are proton-rich; 
some are on the proton drip line, and therefore may be 
proton-halo candidates. The measurements were made 
by aiming these projectiles at a stack of thin Si elements 
which served as both targets and detectors. We also re- 
port <jp for most of the projectiles, the exceptions being 
9 C, whose Op we previously reported Q, and 6 Li and 
10 C, which are unbound following single proton removal. 
We report 02 P for the 2-proton-halo candidate 17 Ne [8|, 
which is also unbound after p-removal. A preliminary 
account of this experiment appeared earlier [jj. Those 
data are superseded by the present measurements since 
we have now more effectively identified and rejected the 
unwanted components of the beams used in the experi- 
ment. 

Proton and heavy ion radiotherapy is now utilized in 
the treatment of cancer [lOj. These beams produce nu- 
clear fragments, including short-lived isotopes such as 
11 C, and these can contribute significant dosage. Experi- 
mental or data are needed for both stable and short-lived 
light ions to better quantify both the dosage delivered by 
these fragments and their spatial distribution along the 



beam path. 

We interpret our data using phenomenological strong 
absorption models, as well as more elaborate models such 
as the optical limit of Glauber multiple-scattering theory. 
Breakup reactions are described in an extended Glauber 
model, incorporating second order noneikonal corrections 
and shell model spectroscopic factors, as well as double 
folding optical potentials generated by complex, density- 
dependent and energy-dependent nucleon-nucleon (NN) 
effective interactions. 

Section II of this paper describes the beams, target- 
detector system, and procedure utilized for the mea- 
surements. Section III explains the selection of reaction 
events and the methods used to determine or and o p . 
The theoretical calculations are presented and compared 
with our data in Section IV. Section V briefly summarizes 
our results and the conclusions we draw. 



II. THE EXPERIMENT 

Our measurements of or are similar to those made 
earlier for the He and Li isotopes on Si 01 . Collimated 
monoenergetic projectiles are selected in a multi-detector 
Si telescope. Reaction events are then identified as those 
whose total energy loss in the telescope differs from that 
of non-reacting projectiles. Since the projectile's energy 
decreases as it travels through the telescope, we obtained 
cross sections at different energies by identifying the de- 
tector in which a reaction occurred; i.e., the first detector 
to give a signal different from projectile signals. 

To measure o p , the heavy fragments from proton re- 
moval reactions were observed in particle identification 
spectra from two detectors - the one which stops the 
fragment and the one which precedes it - following those 
in which the breakup reactions occur |7j. 

An 80 MeV/nucleon 20 Ne beam at the National Su- 
perconducting Cyclotron Laboratory was fragmented on 
a 300 mg/cm 2 9 Be target which produced secondary 
beams. These secondary beams were transported by 
the A1200 analyzing system [12] to the Si detector tele- 
scope. In this system the beams were partially purified by 
polyethylene wedges, and their energy dispersions were 
limited to <1% FWHM by analyzing slits. For exam- 
ple, one A1200 setting delivered usable quantities of 11 C, 
12 N, and 13 at energies of 41, 46, and 51 MeV/nucleon, 
respectively. 

The Si target-detector telescope is shown schemati- 
cally as an inset to Figure 1. The first two detectors, 
which we call Dl and D2, were position-sensitive detec- 
tors (PSD's), with thickness 0.2 mm, lateral dimensions 
2.5 cm x 2.5 cm, and 60 cm separation. They measured 
projectile coordinates normal to the beam axis, defining 
beam of 5 mm radius. Energy losses in these detectors 
identified the desired projectiles, but allowed some con- 
taminants as we discuss later. The next five detectors, 
called D3 through D7, were about 0.5 mm thick and had 
active areas from 300 to 600 mm 2 . They stopped all non- 
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FIG. 1: (color online) Scatter plot of AE2 vs. total energy deposited in the telescope, for incident 48 MeV/nucleon lj O. Most 
contaminants ( 16,lr O, 17 Ne, etc.) are removed by gating on AE2; the treatment of 15 N is described in the text. The inset 
shows the target-detector telescope used for measurements. 



reacting projectiles with A>11 and served as targets for 
reactions, including breakup. The final Li-drifted detec- 
tor, D8, had 5 mm thickness and about 4 cm diameter. 
It stopped projectiles with A<10 and the longer-ranged 
breakup fragments. All detectors except Dl were close- 
packed. Detectors D9 and D10 were not used in this ex- 
periment, but were utilized in a concurrent measurement 
of (Tfl's in the Be isotopes . 

The scatter plot of Fig. fj] illustrates our secondary 
beam composition. It shows AE%, the energy deposited 
in D2, vs. the total energy loss in the telescope, for 
properly-aimed projectiles tentatively identified as 15 O 
by their energy loss in Dl. The non-reacting projectiles 
are the most dense group, labeled 15 O. Events in the 
horizontal tail to the left react downstream of D2, and 
the few to the right are either pileup or positive Q-value 
reactions. Events in the vertical tail are non-reacting 
projectiles with significant energy straggling in D2. 

Several contaminant groups, with magnetic rigidities 
in the A1200 system identical to that of 48 MeV/nucleon 
15 0, also are observed in Fig. fl] Most of these, e.g. 16 - 17 
and 17 ' 20 Ne, are so well resolved that a gate on AE 2 re- 
moves them completely. The treatment of the remaining 
15 N contaminant is described in the next section. 



III. EVENT SELECTION AND DATA 
ANALYSIS 

A. an Measurements 

Spectra of total energy deposited in the Si telescope 
by incident 48 MeV/nucleon 15 projectiles appear in 
Fig. |3 Each of these four spectra shows a prominent 
peak for non-reacting projectiles and, below this peak, a 
continuum due to reactions. 

Events in the spectrum labeled Di_2 are from projec- 
tiles with acceptable positions and normal 15 O energy 
losses in Detectors Dl and D2. The 15 reactions be- 
neath the 15 N contaminant peak are extracted by assum- 
ing a linear energy dependence of the 15 O reaction yield 
near that peak. The contribution from small Q-value 
reactions is included by extrapolating the 15 O reaction 
yield just below the full energy peak to its center. The 
probability 772 of a reaction occurring anywhere down- 
stream of Detector D2 then equals the ratio of reactions 
to total events in this spectrum, which are mainly in the 
sharp peak. The spectrum labeled Di_3 has an addi- 
tional gate requiring normal 15 O ionization in D3, and 
therefore gives the probability 773 for a reaction to occur 
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beyond that detector. Since a non-reacting 15 projectile 
drops in energy from 44 to 38 MeV/nucleon in Detector 
D3, the average cross section for this energy range is 

°R = (m ~ m)/ N 3 (1) 

where Detector D3 has N3 target nuclei per unit area. 
Similarly, the <tr's for the energy ranges 38 to 31 
MeV/nucleon and 31 to 22 MeV/nucleon are determined 
by (f?3 — 774VN4 and (774 — 775) /N5, respectively. Cor- 
rections were made for attenuation of the beam in the 
individual detectors; these were less than 1% in all cases. 

Our <jr data are presented in Table [I] The main 
experimental uncertainties are statistical in origin: the 
counts in the reaction region, subtraction of contaminant 
events, and extrapolation to the non-reacting-projectile 
peak center in Fig. [21 Events above the peak are believed 
to be mainly pileup. If treated similarly to those below 
the peak, they would add about 1.5%, or .02 b, to the 
ctr's quoted in Table [J Given the uncertain origin of 
these events, we make no correction to or but add 2% 
uncertainty in quadrature to the experimental uncertain- 
ties. A further 2% uncertainty in the detector thicknesses 
is also included. 

The contaminants produced very few reactions which 
were counted in the reaction continuum. For example, 
there are about 1700 non-reacting 15 N events in the Di_2 
spectrum of Fig. |3 compared with 3.2xl0 6 15 0's. Hence 
the reactions by contaminants add about 0.05% to the 
measured ctr's. No correction was made for this effect 
since it is so small and the contaminants' reaction cross 
sections are not well known. 



B. Up Measurements 

The heavy fragments from proton removal generally 
have longer ranges than the projectiles which produce 
them. For example, all projectiles with A > 11 stopped in 
Detector D7, but most produced some fragments which 
stopped in D8 along with others which stopped in D7. 
The exception was 13 O which had detectable yield only 
for 12 N stopping in D7. Projectiles with A<10, and all 
of their fragments, stopped in D8. Particle identifica- 
tion (PID) was used to select fragments of interest. The 
PID parameter for fragments stopping in Detector 7 is 
obtained from the energy losses AEq and AE 7 , through 
the equation 

PID = A[(AE 6 + AE 7 ) P -(AE 7 ) P ], (2) 

where p — 1.73 and the arbitrary constant A is chosen 
to place the peaks of interest in convenient channels. A 
similar equation, relating energy losses AE 7 and AEg, 
identifies those fragments stopping in D8. 

Fragment identification was basically restricted to 
those produced in Detectors D3 through D5. A tight 
AE2 gate limited the production region in D2 to about 



its last 0.03 mm, which increased the effective production 
target thickness. Fragments produced earlier in D2 were 
rejected because of their decreased ionization. Produc- 
tion in D6 and beyond was excluded by counting only 
events with significantly lower AE5 than that of non- 
reacting projectiles. This rejected some events produced 
in D5 near its rear face, and we made a model correction 
for this loss. 

As an example, we discuss the boron isotopes produced 
by incident 37 MeV/nucleon 11 C. A PID spectrum for 
fragments stopping in D7 appears in Fig. |3J The domi- 
nant isotope produced is 10 B, with lesser amounts of n B 
(through charge exchange) and 8 B (via some more com- 
plex process). The gap between the 8 B and 10 B groups 
reflects the absence of 9 B, which is unbound. 

Several corrections were made to the observed frag- 
ment yields. The largest of these was for rejection of 
fragments produced in D5, near its rear face, by the 
software gate on AE5 mentioned earlier. This correc- 
tion required the fragment momentum distribution in 
the kinematic center of mass system (cm) of the inci- 
dent projectile, which was calculated from Goldhaber's 
[7L fl4| theory. Typically, 10% of the fragments stopping 
in D7, but none of those stopping in D8, were rejected. 
No fragments produced before D5 were rejected by this 
gate; thus, the overall event loss it caused was typically 
3%. Since fragment production in D2 increased the ef- 
fective target thickness by about 2%, the observed yield 
was reduced by this amount. The loss of projectiles by 
reactions in D3 through D5, added to that of produced 
fragments lost to secondary reactions and hence misiden- 
tified, was about 1%. Possible corrections for fragments 
transmitted through D8 or stopping short of D7 were 
found to be negligible. 

The beam contaminants mentioned earlier had no ef- 
fect on these measurements due to their short range. For 
example, the 37 MeV/nucleon 15 N contaminant in the 
48 MeV/nucleon 15 beam (see Figs. ^ and |2J stops in 
Detector D6. In contrast, the 14 N fragments produced 
from 15 breakup reached at least D7. 

The ratio of the combined corrected fragment yield to 
the incident flux (projectiles with acceptable positions 
and energy losses in Dl and D2) gave the breakup cross 
section; these are presented in TablelTTl The uncertainties 
listed therein were found by adding in quadrature the 
statistical uncertainties and an uncertainty of 1/3 of the 
combined corrections to the data. 



IV. THEORETICAL PREDICTIONS, AND 
COMPARISON WITH MEASUREMENTS 

A. Reaction cross sections 

Experimental reaction cross section data are shown in 
Fig. 21 We first compare these data with the predictions 
of the strong absorption model of Kox [l^ . In this model 
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FIG. 2: Total energy spectra of 48 MeV/nucleon 15 incident upon the Si telescope, gated to exclude events which react 
before detectors D3, D4, D5, or D6. The difference in reaction probabilities between curves labeled Di_2 and Di_3 determines 
<jr at the energies the incident projectiles have in D3, after 15 N contaminant events are subtracted. 



the reaction cross section is 
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-Ec, 



(3) 

where A p ^ are the projectile (target) mass numbers, a = 
1.85 is a mass symmetry parameter related to the volume 
overlap of projectile and target, and C(E) is a correction 
related to the transparency of the optical potential. We 
adopt here the linear approximation of Mittig et al. [Tf^ | , 
C(E) = 0.31 + 0.QIA1E/A which gives reasonable results 
for the range of energies studied here. Bending of trajec- 
tories in the target Coulomb field is taken into account 
by the last factor in Eq.0 There, the Coulomb potential 

is evaluated at R c = 1.07 (Ap /3 + A t 1/3 ) + 2.72 fax 

The Kox formula gives excellent results for stable nu- 
clei when the reduced strong absorption radius is fixed at 
ro=l.l fni, and therefore any significant departure from 
its predictions may disclose a halo structure. The calcula- 
tions shown in Fig. ^were done using the standard value 
for 7"o, except for 6 Li and 9 C where a somewhat larger 
value (ro=1.15 fm) was needed to fit the data. This vari- 
ation must relate to the weak binding of the last two 
nucleons in each of these nuclei. It should be noted that 
the variation Aro=0.05 fm adds 130 mb to an for 6 Li at 



30 MeV/nucleon. 

In Figure|3we display the A-dependence of the Kox an 
predictions, using ro = 1.10 fm for all projectiles, at an 
energy of 35 MeV/nucleon. To compare them with our 
measurements, predictions for all projectiles were indi- 
vidually renormalized to best fit the energy dependence 
of each projectiles's an data. The renormalized calcu- 
lations at 35 MeV/nucleon are plotted as interpolated 
"data" in Figure |SJ The smooth A dependence of the 
predictions reproduces the data on average, but is un- 
able to explain the scatter of the reaction cross sections 
observed in the experiment. This observation suggests 
that a more sophisticated theory is needed to explain 
these data. 

We next use the optical limit (OL) of the Glauber mul- 
tiple scattering theory. The reaction cross section is given 
by 



a R = 2tt I b[l - T(b)]db, 
la 



(4) 



where b is the impact parameter, and the transparency 
function (or elastic survival probability) T(b) is given by 

T(b) = exp[- X (6)]. (5) 

The quantity [1 — T(b)] is called the opacity or nuclear 
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FIG. 3: Particle identification spectra for bound B isotopes produced in Detectors D3 through D5 and identified in D6 and D7, 
with 41 MeV/nucleon n C incident upon the telescope. Note the absence of a group due to 9 B, which is unbound. 



shadow. The scattering phase is given by 




(6) 

where the sum runs over all NN isospin channels. Here 
we assume a zero-range nuclear force and charge symme- 
try for the free NN cross sections (a pp — a nn ). These are 
taken from the parametrization of John et q/.|l7|. We 
ignore in-medium effects, Pauli blocking, and Fermi mo- 
tion, and assume a purely imaginary forward NN scat- 
tering amplitude. In principle, one should correct Eq. 
0] for Coulomb dissociation effects. Estimation of this 
mechanism (see next subsection) leads to the conclusion 
that its contributions are negligibly small. They amount 
to 2 mb for 7 Be and 20 mb for 12 N Coulomb breakup 
at 35 MeV/nucleon, well within the experimental uncer- 
tainties. The profile functions p needed in Eq. [S]are cal- 
culated by Abel transformation of the ordinary particle 
densities p in coordinate space, 



P(b)= / p(Vb 2 + z 2 )dz (7) 

J — oo 

Correction for bending of the trajectories in the target 
Coulomb field is introduced as follows. Eq.0|is rewritten 

as 

/•OO 

<T R = 2vr / b[l - T{b')]db (8) 
Jo 



where b' — (X/k) \ r] + \/t] 2 + k 2 b 2 j is the impact pa- 
rameter for a grazing Coulomb trajectory, ry is the Som- 
merfeld parameter, and k the wave number. After some 
algebra one obtains 

V-pL-Ml). (9) 

Using V c (b') ss V C (R C ) with R c defined above, one ob- 
tains 

a R = 2tt ( l - Y^M) [°° y {1 T{h ' )]db > ( io) 

which also justifies the correction made in Eq. J2J. 

The single particle densities used in this analysis were 
obtained from a standard spherical HF+BCS calculation 
using the density functional of Beiner and Lombard • 
The surface strength of the functional has been slightly 
adjusted to reproduce the known experimental binding 
energies. The rms radii from HF calculations are listed in 
Table ITTT1 They show a surprisingly good agreement with 
experimental data from high energy reactions |l9j | and 
those obtained from the present an data, as described 
below, especially for the loosely bound nuclei 7 Be, 9 C, 
12 N, 13 O and 17 Ne. The results with the OL model are 
displayed in Figs. 0] and 03 Clearly, the mass dependence 
of the reaction cross section is well reproduced since this 
approximation incorporates realistic densities. 
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FIG. 4: (color online) Experimental reaction cross sections (filled points) compared with predictions from the Kox phenomeno- 
logical formula (dash-dotted lines), the optical limit of the Glauber model (dashed lines) and the JLM microscopic folding 
model (continuous lines). Open symbols are data from our earlier work [TlLllflj . See text for details. 



In the remainder of this section we discuss the abil- 
ity of the JLM folding model to describe reaction 
cross sections as a further check of our densities. We 
adopt their nuclear matter approach which incorporates 
a complex, energy- and density-dependent parametriza- 
tion of the effective interaction obtained in the Brueckner 
Hartree-Fock approximation from the Reid hard-core NN 
potential. Studies of elastic scattering of p-shell nuclei 
[2ll l22| indicate that the absorptive component of the 
JLM potential is realistic for loosely bound nuclei and 



needs no renormalization ( N w sal, see below ), while the 
real part needs a significant renormalization. Here we 
extend these studies much closer to the proton drip line. 

In the JLM model the complex form factor for the 
optical potential is given by 

U(R) = J df p drtPp(r p )pt(r t )v(p,E,s) (11) 

where v is the (complex) NN interaction, p p M are the 
single particle densities of the interacting partners, s = 
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FIG. 5: (color online) Interpolated reaction cross sections at 35 MeV/nucleon (points) as a function of projectile mass number. 
Theoretical predictions use the Kox model with the standard reduced strong absorption radius ro = l.l fm (dash-dotted line), 
Glauber (OL) model (dashed line) and JLM model (full line). Open symbols are optical model values from elastic scattering 
of 6 Li+ 28 Si at 318 MeV (triangle, 0) and from 7 Li+ 28 Si at 350 MeV (circle, |2j|). 



r p + R — r t is the NN separation distance between inter- 
acting nucleons and p is the overlap density. The effective 
NN interaction contains an isovector component, which 
gives a negligibly small contribution for p-shell nuclei. 
However, it is included here for convenience in conjunc- 
tion with appropriate single particle isovector densities. 
The coupling of the entrance channel to the breakup and 
particle transfer reactions has been described by a dy- 
namic polarization potential (DPP) which is strong and 
has complicated dependence on radius, mass, and energy 
|23| . To simulate the radial dependence of a DPP, and to 
increase the flexibility of the folding potential we intro- 
duce a smearing function h(r) to obtain our final folding 
potential, 



U(R)= dR'U(R')h(\R- R'\) 



(12) 



The smearin g fu nction h(r) is taken as a normalized 
Gaussian jUl^, 



h(r) 



1 



t 3 n 3/2 



exp(-r^/r) 



(13) 



which behaves as a (5-function for t — ► 0, while for fi- 
nite t values it modifies the rms radius of the folding 
form factor by r\ = (3/2)i 2 , leaving the volume integral 
unchanged. It turns out that the smearing procedure de- 
scribed above is essential in simulating the complicated 
radial dependence of the dynamic polarization potential 
[2lj| . To be consistent with the JLM model we take the 
overlap density in Eq. Illl to be given by 



1 



1 



P P (r p + -zs)Pt{r t - -s) 



1/2 



(14) 



This approximation is physically reasonable since the 
overlap density approaches zero when one of the inter- 
acting nucleons is far from the core, and approaches the 
nuclear matter saturation value for complete overlap. We 
recall that the JLM model was developed to describe the 
optical potential for a nucleoli traversing nuclear matter, 
and its density dependence is defined for densities not 
exceeding the saturation value in nuclear matter. 

The model contains four parameters: two normaliza- 
tion constants N V ,N W and two range parameters t v ,t w . 
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They have been fixed here close to standard values for 
p-shell nuclei [HHI]: ^=0.4, N w =0.85 and t v = t w = 
1.2 fm. We assume that the energy dependence is weak, 
and use this set of parameters at all energies and for all 
projectiles. We are merely interested in a general assesse- 
ment of the JLM model rather than fitting a particular 
cross section. The results of this approach are shown 
with continuous lines in Figs. 0]and|3J Fig. [3] also in- 
cludes the optical model cross sections obtained by fitting- 
known elastic scattering distributions 0, E3 • One ob- 
serves an even better description of the scatter effect in 
the reaction cross sections. 



B. Nuclear radii 

It is interesting to use reaction cross section data to 
infer nuclear sizes. At the energies of our experiment, 
the total NN reaction cross section is large, the mean 
free path is small, and therefore the heavy ion reaction 
cross section is strongly influenced by NN collisions at 
the surface. We chose to describe the tails of the nuclear 
matter densities by a Gaussian [2^: 



Pi{r) = p l0 exp(-r 2 /a 2 ), i=p,t, 



with the normalization, 



PiO 



(15) 



(16) 



The slope parameter is related to the nuclear size by 



1 ^ o 

< rf >= -at 



We let a 1 = a p + af and 



_ 2 a p a t 

Xo — &NNn PpOPtO a 2 ' 



(17) 



(18) 



where ct/vat is the isospin-averaged NN cross section. Un- 
der these assumptions, Eq. (@J can be solved analytically, 



<tr = ita 1 (C + lnxo + £i(xo)) 1 



V b 



(19) 



where C = 0.5772 is the Euler constant and E\ is the 
exponential integral. Eq. 119|) shows essentially the ge- 
ometric character of the reaction cross section since the 
leading term is <7r ~ 7r ( r p + r t)- The leading term in 
the Kox formula is <tr f=a 7r(r p + r f ) 2 . The energy depen- 
dence is governed entirely by a n n- To extract nuclear 
sizes from experimental data, we use Ea.l[T9"|. with the 
target rms radius fixed to r t =3.05 fm which we obtained 
from HF calculations. The compilation of Angeli [27j for 
rms charge radii indicates a value of r c ^=3.12 fm for 28 Si. 
We calculate an upper limit of 2% uncertainty from this 
difference in our theoretical estimations. In practice, the 
parameter a p is gridded in small steps until the calculated 



cross section equals the experimental value. The results 
are listed in Tabic UTTI and displayed in Fig. |SJ These val- 
ues are consistent with those extracted from high energy 
data within their uncertainties. The values and uncer- 
tainties quoted in the tabic include the weighted averages 
over various energies measured here, and statistical un- 
certainties. No provision has been taken to correct values 
in Table III for nucleon finite size. Assuming a nucleon 
rms value of r 2 «0.8 fm 2 , this would increase values in 
Table III by »5%. 

The sensitivity of the cross section to the functional 
form of the single-particle density is studied using har- 
monic oscillator wave functions appropriate for the pd- 
shell nuclei with s and d state admixtures. 



Pr{r) 



1 



— n ( ]X 
15 



exp (- 



where n s 



(20) 

and rid are occupation numbers in the s, 
p and d shells, b T is the range parameter and x = r/b T . 
The finite range of the nuclear force is included in Eq. (JSJ) 
by replacing the S force by a finite range NN interaction 



v(b) = exp (~b 2 1 [i 1 



(21) 



with the interaction range /i=l fm |28|. The normaliza- 
tion in Eq. I|21|l ensures that the reaction volume overlap 
in Eq. EJis not changed, but the smearing enhances the 
weight of the density tail and reduces the contribution 
from the central part of the density. In addition, the 
number of NN inelastic scatterings in the overlap vol- 
ume is weighted differently by a pp a np . As a result, 
the rms radii extracted with this method (shown in Ta- 
ble are slightly smaller than with the Karol model. 
Woods-Saxon (WS) shapes are more appropriate for test- 
ing the role played by long tails in the projectile neutron 
and/or proton densities 



Pr{r) = PQrl 



1 + exp 



p,n 



(22) 



where the po T are normalization constants. Equations 
(1 1171 have been solved numerically. The half-radius pa- 
rameter R T is the most sensitive in extracting the nuclear 
rms radius. No anomalies were found in proton or neu- 
tron surface thickness when using a p m a n ~0.5 fm for 
most nuclei, and slightly larger values for 9 C, 12 N and 
17 Ne. The uncertainties shown in Table ITTT1 are evaluated 
with the bootstrap method of Efron 29] , which has the 
merit of weakening the influence of systematic errors in 
the data, and leads to meaningful ranges in the extracted 
parameters. All these calculations lead to mutually con- 
sistent rms radii for the nuclei studied here ( Table ITTH) . 



C. One- proton removal 

We use the core spectator model and Glauber multiple 
scattering theory to calculate one-proton breakup cross 
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FIG. 6: rms radii extracted from reaction cross sections measured in this experiment (filled symbols) using the gaussian 
density (filled circles), harmonic oscillator density (squares) and Woods-Saxon density (triangles)) compared with the values 
extracted from high-energy interaction cross section [l9l| (open circles). The line shown is 1.075A 1//3 fm to guide the eye. Error 
bars are la deviation arising from energy dependence and statistics. 



sections. This model has been tested extensively on a 
number of one-neutron removal reactions of neutron-rich 
nuclei in the psd shell [30j . We approximate the ground 
state of the projectile (J 7r ) by a superposition of con- 
figurations of the form [I£ ® nlj] J , where I c denotes 
the core states and nlj are the quantum numbers speci- 
fying the single particle wave function of the outermost 
proton. The s.p. wave functions were obtained for a 
Woods-Saxon potential using the effective separation en- 
ergy recipe Spf* = S p + E£ x , where E° x is the excitation 
energy of the core state. The depth of the potential was 
varied in such a way that the known one-proton separa- 
tion energy was reproduced. The radius of the WS po- 
tential was fixed to values close to the rms radius of the 
core (Table ^Qj and the diffuseness was fixed to a = 0.5 
fm for all cases. 

Fine tuning of these parameters could significantly im- 
prove the theoretical results since they determine the 
asymptotic normalization coefficient of the outermost 
proton. We neglect dynamical excitation of core states 
in the reaction. In this approximation, the reaction can 



populate core states only to the extent that there is a 
nonzero spectroscopic factor C 2 S(I c ,nlj) in the projec- 
tile ground state. The various configurations are assumed 
to contribute incoherently to the total breakup cross sec- 
tion for a given core state, 

= E ° 2S ^^ nlj)cr sp (nlj, S e /f). (23) 

nlj 

The total breakup cross section (a p ) is then the sum over 
all particle-bound states of the core. In the present exper- 
iment, the core states were not identified and the removed 
proton was not detected. Therefore, the dynamical fac- 
tor a sp in Eq. (|23|l includes contributions from stripping, 
nuclear dissociation and Coulomb dissociation. Other 
more violent channels, such as core breakup, have been 
ignored. The nuclear mechanisms have been described 
using transition operators defined in terms of scattering 
functions (S-matrix) for p-target and core-target interac- 
tions generated in the JLM folding model as described 
above. Scattering functions in impact parameter repre- 
sentation have been calculated in the eikonal approxima- 
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tion, including noneikonal corrections up to second order. 
Coulomb dissociation has been calculated in the first or- 
der of perturbation theory. Both El and E2 amplitudes 
have been included. The interference of these amplitudes 
does not contribute to the inclusive removal cross section. 

The spectroscopic factors C 2 S(I c ,nlj) employed here 
were calculated with the shell model code OXBASH [2l| 
using the WBP 32] interaction within a Is — lp— 2slrf — 
l/2p model space. Excitations up to 4hu) have been in- 
cluded in this single particle model space for most of the 
cases. Inclusion of higher nhuj excitations with this effec- 
tive interaction was shown to better describe some exper- 
imental observables |3^| . Spurious center-of-mass compo- 
nents of the projectile have been suppressed by the usual 
method (sij of adding a center-of-mass hamiltonian to 
the nuclear interaction. In all cases, the experimentally 
established spin-parity and core excitation energies have 
been used. The shell model spectroscopic factors have 
been multiplied by the center-of-mass correction factor 
A P /(A P -1), following Refs. [M |H . 

The projectile ground state spin-parity (J 7r ), core spin- 
parity and excitation energies (E ea; , LT), the stripping, 
diffractive dissociation and Coulomb dissociation cross 
section, as well as the shell model spectroscopic fac- 
tors are listed in Table IIVI The total inclusive cross 
section a Giauber wag correc ted for center-of-mass effects 
as explained above. All calculations were done at 35 
MeV/nucleon. 

Comparison with experimental values in TablelTTlshows 
reasonable agreement. It should be noted that the ex- 
tended Glauber model is designed to describe one nucleon 
removal reactions for loosely bound nuclei. The wave 
function of the outermost nucleon is assumed to pene- 
trate substantially into the classically forbidden region, 
since most of the reaction probabilty is localized at the 
surface. The interaction should be strongly absorptive. 
In addition, the core is assumed to be compact and well 
decoupled from the outermost nucleon. A special case 
is 10 B which has a quite large one-proton separation en- 
ergy (S p 6.8 MeV), and a relatively fragile core, 9 Be. In 
this case the reaction is not completely peripheral, with 
significant contributions arising from impact parameters 
b < R c , where R c is the core radius. The calculated 
breakup cross section exceeds the experimental value by 
more than 50%. Only by excluding contributions from 
the nuclear interior can one obtain a reasonable value 
(Table El. 



D. Discussion of Results 

The reaction cross sections for 6 Li measured in this 
experiment and at higher energy |llj | show in Fig. 01 a 
smooth energy dependence, falling at 50 MeV/nucleon, 
which is well fitted by the Kox model and slightly un- 
derpredicted by the Glauber models. For 7 Be we see a 
flat energy dependence, with all models slightly under- 
predicting the 40 MeV/nucleon datum. 



Figures 0] and [3] show a large reaction cross section for 
the 9 C projectile. Both Glauber (OL) and JLM models 
predict this increase, in reasonable agreement with the 
experiment. The rms radius obtained from our present 
data exceeds by about 12% the value obtained from high 
energy data [lj|, and it is close to that obtained for 17 Ne, 
another loosely bound nucleus. The smaller Coulomb 
barrier in 9 C enhances the penetrability of the outermost 
proton much beyond the range of nuclear forces and this 
compensates the larger interaction volume provided by 
17 Ne. Unfortunately, the uncertainty is relatively large 
and the halo structure of this nucleus cannot be firmly 
established. More accurate experimental data are needed 
to clarify this point . 

The reaction cross sections for 10 B, 10 C, and n C are 
definitely smaller than for 9 C and not much larger than 
those of 6 Li and 7 Be. The extracted rms radius of n C is 
quite small though in agreement with high energy data 
|19| . An independent measurement by Liatard et q?.|37j 
indicates a much larger value, 2.46±0.30 fm, in agreement 
with our HF calculation of 2.48 fm (see Table EH). 12 N 
seems to behave normally. Both Glauber (OL) and JLM 
models give an excellent description of the cross section. 
The extracted radii seem to be in good agreement with 
HF calculations (see Table llTl|) . 

Figure [5] shows that or stays nearly constant from 12 N 
through 15 0, as predicted by the JLM and OL (but not 
the Kox) models. For well-bound projectiles we expect 
an increase with increasing A. However, S p is 0.6, 1.5 
and 7.3 MeV for 12 N, 13 and 15 0, respectively, and 
the weaker binding of the two lighter nuclei increases the 
range of their valence protons. The two effects appear to 
roughly compensate for each other. 

Among the nuclei we studied, 12 N has by far the largest 
one proton removal reaction cross section. Furthermore, 
the p-removal process for this case is related to 11 C+p 
radiative capture. Depending on their initial CNO abun- 
dances, this reaction may have been important in some 
super-massive stars in the early universe, allowing the 
stars to explode as supernovae rather than collapsing 
as black holes before ejecting any mass.|38l |3JJ. The 
asymptotic normalization coefficient (ANC) for 12 N — > 
11 C+p has been measured recently using transfer reac- 
tions [40j. However, these reaction-rate data are ambigu- 
ous and therefore it is desirable to remeasure the ANC 
by breakup, as an independent check. The one proton 
separation energy is small (S p — 0.6 MeV), and thus 12 N 
could be a proton halo candidate. The Glauber model 
nicely reproduces the measured a p . However, the shell 
model calculations suggest a very fragmented structure 
for the g.s. wave function, significant contributions aris- 
ing from (g.s, 3/2-®lp 1/2 ), (g.s, 3/2"(g)lp 3/2 ) ^=2.0 
MeV, l/2-(g)lp 3/2 ) and (E x =4.8 MeV, 3/2-®lp 1/2 ). 
The ANC's have been measured |4£|, for the core ground 
state components; they are C 2 Pl/2 = 1.4 ± 0.2 fm -1 and 
C 2 p = 0.33 ± 0.05 fm -1 . Since the single particle nor- 
malization coefficients are almost identical for the lpi/2 
and IP3/2 wave functions, one should have approxima- 
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tively S lpi/2 /S lp3/2 « C? Pl/a /C? Pa /a' a relation that is 
satisfied for the states mentioned above for which the ex- 
perimental values of the ANC are known, when one uses 
the corresponding shell model spectroscopic factors listed 
in Table IIVI The reaction cross section measurements 
showed no anomaly for this nucleus. The calculated par- 
allel and transverse core fragment momentum distribu- 
tions are narrow (FWHM 2 - 86 MeV/c, FWHM, - 120 
MeV/c) comparable to that found for a well established 
halo nucleus, such as 8 B |4|. The halo character of 12 N 
could not be firmly established on the basis of the present 
data. More precise reaction and breakup cross section 
data are needed as well as a clear separation of contribu- 
tions from core excited states. 

For halo nuclei, one can assume that the core is decou- 
pled from the halo nucleon, and the following decompo- 
sition of the absorption operator (Eq. 0} holds 

1 - T « 1 - T c T fc = (1 - T c ) + T c (l - T fc ) (24) 

Here, the first bracketed term on the right describes the 
core absorption (in the absence of the halo particle), 
while the second describes the absorption of the halo par- 
ticle, weighted by the core survival probability. One can 
associate this last term with the stripping component of 
the total breakup cross section a p . Assuming that the 
diffractive dissociation and Coulomb dissociation com- 
ponents are small, one can approximate <tr = a R + a p . 
The difference of the ctr's given for n C and 12 N in the 
middle energy bins of Table [Q (at energies 28 and 33 
MeV/nucleon, respectively) is 110 ± 86 mb. This agrees 
with the breakup cross section of 120 ± 6 mb given in 
Table ITU though with a large uncertainty. 

15 O has a one-proton separation energy S p four times 
larger than 13 0, yet the 15 breakup cross section is two 
times larger. We expect a larger S p to give a smaller 
a p and so the Glauber model, which is very sensitive 
to the separation energy, is not able to reproduce this 
behaviour. This anomaly probably results from compe- 
tition with 2p-breakup; we found that the 13 breakup 
produces more C than 12 N, though the C yield is isotopi- 
cally unresolved. Other similar cases were observed for 
nuclei which have one nucleon outside a fragile core. For 
eaxmple, 12 F3e, which has one neutron loosely bound to 
a neutron halo nucleus 11 Be, has a smaller Cnthan ai n 
[l3| . 9 C behaves similarly for proton removal 0). 

V. CONCLUSIONS 

We have measured reaction cross sections ur and one- 
proton removal cross sections a p for a range of stable and 
short-lived nuclei close to the proton drip line. These ur 
extend our earlier measurements to projectiles of higher 
A and Z. Comparison with earlier, higher energy data for 
6 Li and 7 Be shows a reasonable energy dependence. The 
optical limit of the multiple scattering Glauber model, 
and double folding optical potentials derived from the 
NN effective interaction in nuclear matter, are about 



equally successful in describing the data, including their 
energy- and mass-dependence. However, more accurate 
data would be needed to distinguish between them or to 
test sensitively for nuclei with weak halos. 

Root-mean-square nuclear radii have been extracted 
using Glauber theory and three different functional forms 
for the projectile single particle densitities. These results 
are mutually consistent and in good agreement with val- 
ues extracted from high energy data, except for 9 C in 
which case a significantly larger radius has been found. 

One proton removal cross sections vary widely These 
reactions were described within an extended Glauber 
model by incorporating fundamental NN interactions and 
spectroscopic factors from the shell model. The large 
a p in the 12 N case suggests a possible weak halo, al- 
though the wave function is fragmented into many com- 
ponents. An anomaly has been found in the A=12-15 
region, where the one-proton removal cross sections and 
total reaction cross sections do not behave as expected 
from proton separation-energy systematics. More ex- 
perimental and theoretical effort shoud be devoted to 
a better understanding of the competition between — lp 
and — 2p channels for nuclei, such as 13 O, where the last 
proton's separation energy exceeds that of the next one 
(i.e.,S p > 2^2p)- 
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VII. TABLES 



TABLE II: Proton-removal cross sections for light nuclei on 
Si, averaged between the specified energies. The Glauber the- 
oretical predictions are described in Section fl V CI 



TABLE I: Total reaction cross sections, o\r, for ten light iso- 
topes on Si, averaged between the listed energies. 



System 


Energy interval 


OR 




(MeV/nucleon) 


(b) 


Ti+Si 


15.1-19.4 


1.55L0.06 




19.4-23.1 


1.54L0.06 


'Be+Si 


21.0-26.2 


1.54L0.09 




26.2-30.6 


1.44L0.10 




30.6-34.7 


1.54L0.10 


U, B+Si 


25.0-30.0 


1.58L0.06 




30.0-34.3 


1.65L0.07 




34.3-38.4 


1.58L0.07 


9 C+Si 


27.7-34.9 


1.75L0.18 




34.9-41.0 


1.65L0.19 




41.0-46.6 


1.69L0.20 


iu C+Si 


22.6-30.0 


1.53L0.06 




30.0-36.2 


1.69±0.07 




36.2-41.7 


1.57L0.07 




41.7-46.4 


1.52L0.07 


"C+Si 


17.3-25.3 


1.61L0.05 




25.3-31.7 


1.64±0.05 




31.7-37.0 


1.54L0.05 


i;2 N+Si 


20.3-29.2 


1.74L0.07 




29.2-36.4 


1.75L0.07 




36.4-42.2 


1.68L0.08 


ia O+Si 


22.5-32.3 


1.74L0.11 




32.3-40.4 


1.77L0.12 




40.4-46.8 


1.69L0.13 


i;> 0+Si 


22.1-30.8 


1.77L0.06 




30.8-38.1 


1.81L0.06 




38.1-44.0 


1.62L0.06 


iY Ne+Si 


27.5-37.7 


1.99L0.09 




37.7-46.3 


1.92±0.10 




46.3-53.3 


1.84L0.10 



Projectile, 


Fragment 


Energy interval 


Measurement 


Glauber 


cross section 




(MeV/nucleon) 


(mb) 


(mb) 


'Be, CTi p 


Ti 


26-38 


90±6 


117 


£>, &i P 


9 Be 


30-42 


41±3 


62 


iir 1 „ 

L/, (7i p 


io B 


17-37 


53±2 


78 


12 N, <T lp 


U G 


20-42 


120±6 


130 


13 0, a lp 


12 N 


23-49 


31±5 


56 


15 0, a lp 


14 N 


22-44 


64±3 


46 


17 Ne,cr 2 p 


15q 


28-53 


223L18 





TABLE III: Proton (r^), neutron (r„), and matter (r m ) rms 
radii from HF calculations are given in Columns 2 through 
4. Column 5 gives experimental values extracted from high- 
energy reaction cross sections [Til . The last three columns 
show experimental values extracted from our present data, 
using Glauber theory (OL) and assuming gaussian, harmonic 
oscillator, and Woods-Saxon densities. 



Nucleus 




T u 


'771 


r m (exp) [19] 


r m (gauss) 


r m (HO) 


r m (WS) 


b Li 


2.33 


2.31 


2.32 


2.32± 0.03 


2.28± 0.10 


2.23±0.11 


2.30±0.12 


7 Be 


2.49 


2.24 


2.39 


2.31± 0.02 


2.28± 0.20 


2.20±0.19 


2.28±0.20 


10 B 


2.46 


2.46 


2.46 


2.20± 0.06 


2.29± 0.13 


2.24±0.13 


2.27±0.13 


9 C 


2.76 


2.25 


2.60 


2.42± 0.03 


2.75± 0.34 


2.71±0.26 


2.69±0.25 


io c 


2.57 


2.31 


2.47 


2.27± 0.03 


2.32± 0.14 


2.29±0.14 


2.31±0.14 


n C 


2.53 


2.41 


2.48 


2.12± 0.06 


2.20± 0.10 


2.18±0.11 


2.23±0.10 


12 N 


2.70 


2.45 


2.60 


2.47± 0.07 


2.51± 0.13 


2.49±0.12 


2.60±0.12 


HQ 


2.81 


2.47 


2.69 


2.53± 0.05 


2.54± 0.21 


2.54±0.20 


2.59±0.20 


15 


2.69 


2.58 


2.64 


2.44± 0.04 


2.40± 0.10 


2.40±0.10 


2.45±0.10 


17 Ne 


2.98 


2.61 


2.83 


2.75± 0.07 


2.82± 0.15 


2.84±0.14 


2.89±0.14 
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TABLE IV: Calculated spectroscopic factors (C 2 S) and cross 
sections <r(IJ) to the core excited states (E e x,Ic) popu- 
lated in single proton removal from the projectile ( A Z,J' r ) 
by the silicon target. The contribution arising from strip- 
ping a abs , diffractive dissociation (adiff) and Coulomb 
dissociation(cr COU !) are detailed. The total inclusive cross 
section a ^ lauber ; s corrected for center-of-mass effects as ex- 
plained in the main text. 



A Z J* 




I? 


nlj 


C 2 S 


<7 a bs 
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